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ABSTRACT 

In this second part we present a set of methods, analytical and numerical, which can describe behaviour 
in (non) equilibrium ensembles, both classical and quantum, especially in the complex systems, where 
the standard approaches cannot be applied. The key points demonstrating advantages of this approach 
are: (i) effects of localization of possible quantum states; (ii) effects of non-perturbative multiscales which 
cannot be calculated by means of perturbation approaches; (iii) effects of formation of complex/collective 
quantum patterns from localized modes and classification and possible control of the full zoo of quan- 
tum states, including (meta) stable localized patterns (waveletons) . We demonstrate the appearance 
of nontrivial localized (meta) stable states/patterns in a number of collective models covered by the 
(quantum) /(master) hierarchy of Wigner-von Neumann-Moyal-Lindblad equations, which are the result 
of "wignerization" procedure (Weyl-Wigner-Moyal quantization) of classical BBGKY kinetic hierarchy, 
and present the explicit constructions for exact analytical/numerical computations. Our fast and efficient 
approach is based on variational and multiresolution representations in the bases of polynomial tensor 
algebras of generalized localized states (fast convergent variational- wavelet representation) . We construct 
the representations for hierarchy/ algebra of observables(symbols)/distribution functions via the complete 
multiscale decompositions, which allow to consider the polynomial and rational type of nonlincaritics. 
The solutions are represented via the exact decomposition in nonlinear high-localized eigenmodes, which 
correspond to the full multiresolution expansion in all underlying hidden time/space or phase space 
scales. In contrast with different approaches we do not use perturbation technique or linearization proce- 
dures. Numerical modeling shows the creation of different internal structures from localized modes, which 
are related to the localized (meta) stable patterns (waveletons), entangled ensembles (with subsequent 
decoherence) and/or chaotic-like type of behaviour. 
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ABSTRACT 

In this second part we present a set of methods, analytical and numerical, which can describe behaviour in 
(non) equilibrium ensembles, both classical and quantum, especially in the complex systems, where the standard 
approaches cannot be applied. The key points demonstrating advantages of this approach are: (i) effects of 
localization of possible quantum states; (ii) effects of non-perturbative multiscales which cannot be calculated by 
means of perturbation approaches; (iii) effects of formation of complex/collective quantum patterns from localized 
modes and classification and possible control of the full zoo of quantum states, including (meta) stable localized 
patterns (waveletons). We demonstrate the appearance of nontrivial localized (meta) stable states/patterns in 
a number of collective models covered by the (quantum)/ (master) hierarchy of Wigner-von Neumann-Moyal- 
Lindblad equations, which are the result of "wignerization" procedure (Weyl-Wigner-Moyal quantization) of 
classical BBGKY kinetic hierarchy, and present the explicit constructions for exact analytical/numerical compu- 
tations. Our fast and efficient approach is based on variational and multiresolution representations in the bases 
of polynomial tensor algebras of generalized localized states (fast convergent variational- wavelet representation) . 
We construct the representations for hierarchy/ algebra of observables (symbols) /distribution functions via the 
complete multiscale decompositions, which allow to consider the polynomial and rational type of nonlincarities. 
The solutions are represented via the exact decomposition in nonlinear high-localized eigenmodes, which corre- 
spond to the full multiresolution expansion in all underlying hidden time/space or phase space scales. In contrast 
with different approaches we do not use perturbation technique or linearization procedures. Numerical modeling 
shows the creation of different internal structures from localized modes, which are related to the localized (meta) 
stable patterns (waveletons), entangled ensembles (with subsequent decoherence) and/or chaotic- like type of 
behaviour 

Keywords: localization, pattern formation, multiscales, multiresolution, waveletons, (non) equilibrium ensem- 
bles 

1. INTRODUCTION: CLASSICAL AND QUANTUM ENSEMBLES 

In this paper we consider the applications of a numerical-analytical technique based on local nonlinear harmonic 
analysis to the description of quantum ensembles. The corresponding class of individual Hamiltonians has the 
form 

H(p,q) = ^- + U(p,q), (1) 

where U{p,q) is an arbitrary polynomial function on p, q, and plays the key role in many areas of physics [1]. 
Many cases, related to some physics models, are considered in [2]- [8]. It is a continuation of our more qualitative 
approach considered in part I [9] . In this part our goals are some attempt of classification and the constructions 
of explicit numerical-analytical representations for the existing quantum states in the class of models described 
(non) equilibrium ensembles and related collective models. There is a hope on the understanding of relation 
between the structure of initial Hamiltonians and the possible types of quantum states and the qualitative type 
of their behaviour. Inside the full spectrum there are at least three possibilities which are the most important from 
our point of view: localized states, chaotic-like or/and entangled patterns, localized (stable) patterns (qualitative 

http://www.ipme.ru/zeitlin.html, http://www.ipme.nw.ru/zeitlin.html E-mail: zeitlin@math.ipme.ru, an- 
ton@math.ipme.ru 
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definitions/descriptions can be found in part I [9]). All such states are interesting in the different areas of physics 
[1]. Our starting point is the general point of view of a deformation quantization approach at least on the 
naive Moyal/Weyl/Wigner level. The main point of such approach is based on ideas from [1], which allow to 
consider the algebras of quantum observables as the deformations of commutative algebras of classical observables 
(functions). So, if we have the classical counterpart of Hamiltonian (1) as a model for classical dynamics and 
the Poisson manifold M (or symplectic manifold or Lie coalgebra, etc) as the corresponding phase space, then 
for quantum calculations we need first of all to find an associative (but non-commutative) star product * on the 
space of formal power series in h with coefficients in the space of smooth functions on M such that 

f*g = fg + Hf,g} + E s). ( 2 ) 

n>2 

where {/, g} is the Poisson brackets, B n are bidifferential operators. In the naive calculations we may use the 
simple formal rule: 

* = CX P i^C^Q^P - *dp3q)) ( 3 ) 



In this paper we consider the calculations of the Wigner functions W(p,q,t) (WF) corresponding to the 
classical polynomial Hamiltonian H(p,q,t) as the solution of the Wigner- von Neumann equation [1]: 



d 

iti—W 
dt 



H *W -W * H 



(4) 



and related Wigner-like equations for different ensembles. According to the Weyl transform, a quantum state 
(wave function or density operator p) corresponds to the Wigner function, which is the analogue in some sense 
of classical phase-space distribution [1]. We consider the following operator form of differential equations for 
time-dependent WF, W — W(p,q,t): 



W t = — sin 



■HW 



which is a result of the Weyl transform or "wignerization" of von Neumann equation for density matrix: 



(5) 



(6) 



1.1. BBGKY Ensembles 

We start from set-up for kinetic BBGKY hierarchy (as c-countcrpart of proper q- hierarchy) . We present the 
explicit analytical construction for solutions of both hierarchies of equations, which is based on tensor algebra 
extensions of multircsolution representation and variational formulation. We give explicit representation for 
hierarchy of n-particle reduced distribution functions in the base of high- localized generalized coherent (regarding 
underlying generic symmetry (affine group in the simplest case)) states given by polynomial tensor algebra of our 
basis functions (wavelet families), which takes into account contributions from all underlying hidden multiscales 
from the coarsest scale of resolution to the finest one to provide full information about stochastic dynamical 
process. The difference between classical and quantum case is concetrated in the structure of the set of operators 
included in the set-up and, surely, depends on the method of quantization. But, in the naive Wigner-Weyl 
approach for quantum case the symbols of operators play the same role as usual functions in classical case. In 
some sense, our approach resembles Bogolyubov's one and related approaches but we don't use any perturbation 
technique (like virial expansion) or linearization procedures. Most important, that numerical modeling in both 
cases shows the creation of different internal (coherent) structures from localized modes, which are related to 
stable (equilibrium) or unstable type of behaviour and corresponding pattern (waveletons) formation. 

Let M be the phase space of ensemble of N particles (dimM = 6iV) with coordinates 
x i = {q i ,p i ), i = l,...,N, q i = {ql,q 2 i ,ql)&B i , Pi = (ph&Pi) € R 3 ', Q = («i, ■ ■ ■ , Qn) G R 3N ■ (7) 
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Individual and collective measures are: 

N 

Hi = dxi = dqidpi, h = ^ (8) 

i=l 

Distribution function Dn(xi, . . . , Xn', t) satisfies Liouville equation of motion for ensemble with Hamiltonian Hn 
and normalization constraint: 

{H Nl D N }, J D N (xi, . . . 1 x N ;t)dn = 1 (9) 



dt 

where Poisson brackets are: 

N 



Sn n \ \-{dH N dD N dH N dD N \ 

{h n , ^}-L(-^-^- -fa-az-) (10) 

Our constructions can be applied to the following general Hamiltonians: 

N 2 

^ = E(S7 + C/ '(«))+ E Vain,*) (11) 

i=l l<i<j<N 

where potentials Ui(q) = Ui(q±, . . . ,?jv) and Uij(qi,qj) are not more than rational functions on coordinates. Let 
L s and Ljj be the Liouvillean operators (vector fields) 

t - V /"^ fl du i 9 \ V Z Z = ^A + ^)A nol 

s 0« a^J * J % d P i dq dp, y ' 

For s=N we have the following representation for Liouvillean vector field Ln = {Hn, •} and the corresponding 
ensemble equation of motion: 

^L + L N D N = (13) 
Ln is self-adjoint operator regarding standard pairing on the set of phase space functions. Let 

F N (xi, ...,x N ;t) = y^D N (x-i, ...,x N ;t) (14) 

Sn 

be the N-particle distribution function (Sn is permutation group of N element s). Then we have the hierarchy 
of reduced distribution functions (V s is the corresponding normalized volume factor) 

F s (x 1 ,...,x 8 ;t) =V S / D N (x 1 ,...,x N ;t) J| ^ (15) 

J s+l<i<N 

After standard manipulations we arrived to c-BBGKY hierarchy: 

OF 1 f s ^ 

+ L S F S = - d^ s+ i ^ L i,s+iF s+ i (16) 
at v J . =1 

It should be noted that we may apply our approach even to more general formulation. As in the general as in 
particular situations (cut-off, e.g.) we are interested in the cases when 

k 

F k ( Xl , ...,x k ;t) = Y[F 1 (x l ;t) + G k {x 1 , . . .,x k ;t), (17) 

i=l 

where G k are correlators, really have additional reductions as in the simplest case of one-particle truncation 
(Vlasov/Boltzmann-likc systems). So, the proper dynamical formulation is reduced to the (infinite) set of equa- 
tions for correlators/partition functions. Then by using physical motivated reductions or/and during the corre- 
sponding cut-off procedure we obtain, instead of linear and pseudodifferential (in general case) equations, their 
finite-dimensional but nonlinear approximations with the polynomial type of nonlinearities (more exactly, mul- 
tilinearities) . Our key point in the following consideration is the proper generalization of naive perturbative 
multiscale Bogolyubov's structure restricted by the set of additional physical hypotheses. 
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1.2. Quantum ensembles 

Let us start from the second quantized representation for an algebra of observables 



A = (A ,A 1 ,...,A S ,...) 



in the standard form 



A = A a + J dx l ^ + {xi)A 1 ^{x 1 ) + ... + (si) 1 J dx 1 ... dx s ^+(x 1 ) . . . ^+(x s )A s ^(x s ) . . . *(xx) + . . 
N-particle Wigner functions 

f s 1 

W s (xi, . . . ,x s ) = / dkx . . . dk s evp( - iy^^kiPi)Trp^ + (qi - -^hki) ... 

J i=i 

^ + {q s - ^hk s )^(q s + ^hk s ) . . . *(<7i + ^hk 8 ) 

allow us to consider them as some quasiprobabilities and provide useful bridge between c- and q-cases: 
< A>= r TrpA = ^2(s\y 1 / YldniAsixu . . . ,x s )W s {x 1 , . . . ,x s ) 



s=0 



The full description for quantum ensemble can be done by the whole hierarchy of functions (symbols): 

W = {W a (xi,...,x a ), s = 0,l,2...} 
So, we may consider the following q-hierarchy as the result of "wignerization" procedure for c-BBGKY 



d t W s (t,x 1 ,...,x s ) = ^2L° j W s (x 1 ,...,x s ) + ^2^2L j!n W s (x 1 ,...,x s ) 

j—l j<nn= 1 



+ ^2 dx s+ iS(k s+1 )L jtS+1 W s+ i(x 1 ,...,x s+1 ) 

3=1 J 



where 



3 y m' 



dm 



exp {-l^wr^'^^wr^ 



In quantum statistics the ensemble properties are described by the density operator 

p(t) = 2«; j |*i(t) ><*<(t)|, J2 Wi = 1 



After Wcyl transform we have the following decomposition via partial Wigner functions Wi (p, q, t) for the whole 
ensemble Wigner function: 

W{p,q,t) = Y / w i W i {p,q,t) (27) 

i 

where the partial Wigner functions 

W n (q,p,t) = ^ J d^exp(-lpe)*;(9-^,t)*n(«+^,t) (28) 
are solutions of proper Wigner equations: 

ow n P dw n A (-iY(h/2r d^u n ( q ) d 2 ^w n 

8t m dq {2l+l)\ dq 2 ^ 1 dp 2 ^ 1 1 ' 

Our approach, presented below, in some sense has allusion on the analysis of the following standard simple 
model considered in [1]. Let us consider model of interaction of nonresonant atom with quantized electromagnetic 
field: 

v 2 

H = -^- + U(x), U(x) = U Q (z,t)g(x)a+a (30) 

where potential U depends on creation/annihilation operators and some polynomial on x operator function (or 
approximation) g(x). It is possible to solve Schroedinger equation 

dl^ > 

ih^-— = H\V> (31) 
at 

by the simple ansatz 

oo „ 

>= ^2 Wn / Ax ^n{x,t)\x > ®|n > (32) 

-oo J 

which leads to the hierarchy of analogous equations with potentials created by n-particle Fock subspaces 

ih d * n ^ t] = {|| + U (t)g(x)n}* n (x,t) (33) 

where ^ n {x, t) is the probability amplitude of finding the atom at the time t at the position x and the field in 
the n Fock state. Instead of this, we may apply the Wigner approach starting with proper full density matrix 

p=|*(t) ><*(t)|= (34) 



>< x"\ ® In' >< n"\ 



^ w n ,w* n „ f dx' f &x"^ n ,{x',t)^ n „{x",i)\x' 

Standard reduction gives pure atomic density matrix 

Pa= < n\p\n > = (35) 

Jn=Q 

Y,\ w n\ 2 J dx' J ax"^ n {x' ,t)^* n {x" ,t)\x' X x"\ 
Then we have incoherent superposition 

oo 

W{x,p,t) = Y,\wn\ 2 W n (x,p,t) (36) 



n=0 
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of the atomic Wigncr functions (28) corresponding to the atom motion in the potential U n (x) (which is not more 
than polynomial in x) generated by n- level Fock state. They are solutions of proper Wigner equations (29). 

The next case describes the important decoherence process. Let us have collective and environment subsys- 
tems with their own Hilbert spaces 

H = H c ®He (37) 
Relevant dynamics is described by three parts including interaction 

H = H c ®I e +I c ®H e + H int (38) 

For analysis, we can choose Lindblad master equation [1] 

p = 1[JT, p] - ln(L+L nP + P L+L n - 2L nP L+) (39) 

n 

which preserves the positivity of density matrix and it is Markovian but it is not general form of exact master 
equation. Other choice is Wigner transform of master equation and it is more preferable for us 

W = {H, W} PB + 2 2n ( 2 n + 1 )i d " n+1U ^ d p n+lw ^P) + 2 -Td P pW + Dd 2 p W (40) 

n>l ^ 71 >' 

In the next section we consider the variational-wavelet approach for the solution of all these Wigner-like 
equations (4), (5), (6), (23), (29), (40) for the case of an arbitrary polynomial U(q,p), which corresponds to a 
finite number of terms in the series expansion in (5), (29), (40) or to proper finite order of h. Analogous approach 
can be applied to classical counterpart (16) also. Our approach is based on the extension of our variational- 
wavelet approach [2]- [8]. Wavelet analysis is some set of mathematical methods, which gives the possibility to 
take into account high- localized states, control convergence of any type of expansions and gives maximum sparse 
forms for the general type of operators in such localized bases. These bases are the natural generalization of 
standard coherent, squeezed, thermal squeezed states [1], which correspond to quadratical systems (pure linear 
dynamics) with Gaussian Wigner functions. The representations of underlying symmetry group (affine group 
in the simplest case) on the proper functional space of states generate the exact multiscale expansions which 
allow to control contributions to the final result from each scale of resolution from the whole underlying infinite 
scale of spaces. Numerical calculations according to methods of part I [9] explicitly demonstrate the quantum 
interference of generalized localized states, pattern formation from localized eigenmodes and the appearance of 
(stable) localized patterns (waveletons) . 

2. VARIATIONAL MULTIRESOLUTION REPRESENTATION 

2.1. Multiscale Decomposition for Space of States: Functional Realization and Metric 
Structure 

We obtain our multiscale/multiresolution representations for solutions of Wigner-like equations via a variational- 
wavelet approach. We represent the solutions as decomposition into localized eigenmodes (regarding action of 
affine group, i.e. hidden symmetry of the underlying functional space of states) related to the hidden underlying 
set of scales [10]: 

oo 

W n (t,q, P ) = ®W n (t,q,p), (41) 

i— i c 

where value i c corresponds to the coarsest level of resolution c or to the internal scale with the number c in the 
full multiresolution decomposition of the underlying functional space (L 2 , e.g.) corresponding to the problem 
under consideration: 

V c C V c+1 C V c+2 C . . . (42) 
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and p = (pi,p2, ■ ■ .), q = (qi,q 2} ■ • •)' x i = (Pi) 9i) •••iPu Qi) arc coordinates in phase space. In the following 
we may consider as fixed as variable numbers of particles. 

We introduce the Fock-likc space structure (in addition to the standard one, if we consider second-quantized 
case) on the whole space of internal hidden scales. 

# = 00^" ( 43 ) 

i n 

for the set of n-partial Wigner functions (states): 

W* = {Wl Wi(x V ,t), . . . , Wi,( Xl , r v :/:... .}, (44) 

where W p (x\, . . . , x p ; t) G H p , H° = C, H p = L 2 (R 6p ) (or any different proper functional space), with the 
natural Fock space like norm: 

{W^^Wi + Y, [w?(x 1 ,...,x i ;t)f[ fM . (45) 

First of all, we consider W = W(t) as a function of time only W G L 2 (R), via multircsolution decomposition 
which naturally and efficiently introduces the infinite sequence of the underlying hidden scales [10]. We have the 
contribution to the final result from each scale of resolution from the whole infinite scale of spaces (16). The 
closed subspace Vj(j G Z) corresponds to the level j of resolution, or to the scale j and satisfies the following 
properties: let Dj be the orthonormal complement of Vj with respect to Vj + \\ Vj + \ = Vj 0-Dj- Then we have 
the following decomposition: 



{W(t)}= D j= V c ($D j: (46) 

— oo<j<oo j—0 

in case when V c is the coarsest scale of resolution. The subgroup of translations generates a basis for the fixed 
scale number: span fegZ {2- J / 2 4'(2 : 't — k)} = Dj. The whole basis is generated by action of the full affinc group: 

spzn keZjeZ {2 j/2 *(2 J t - k)} = span fcJeZ {* ji k } = {W{t)} (47) 

2.2. Tensor Product Structure 

Let sequence 

{F/}, V; G L 2 (R) (48) 
correspond to multiresolution analysis on time axis and 

{Vp}, VpcL 2 (R) (49) 
correspond to multiresolution analysis for coordinate x i7 then 

Vf +1 = Vf 1 ® . . . Vf" Vf (50) 
corresponds to multiresolution analysis for n-particle distribution fuction W n (xi, . . . , x n ;t). E.g., for n = 2: 

v o = {/ : f( x i, x 2) = a kite4> 2 {xi - h,x 2 - k 2 ), a kuk2 G £ 2 (Z 2 )}, (51) 

ki ,k 2 

where 

<j> 2 (x 1 ,x 2 ) = 4> 1 {x 1 )4> 2 {x 2 ) = (j) 1 ® 4> 2 { Xl ,x 2 ), (52) 
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and 4> l (xi) = <j){xi) form a multiresolution basis corresponding to {Vp}. If 

{^(si -£)}, lEZ (53) 

form an orthonormal set, then 

4> 2 (x 1 -k 1 ,x 2 -k 2 ) (54) 

form an orthonormal basis for Vq. Action of affine group provides us by multiresolution representation of L 2 (R 2 ). 
After introducing detail spaces D 2 , we have, e.g. 

V 2 ^V 2 ®D 2 . (55) 

Then 3-component basis for D 2 is generated by translations of three functions 

^ 2 1 = <j> 1 (x 1 )^>^ 2 (x 2 ), 

^ 2 = ^ 1 (x 1 )^<p 2 (x 2 ), (56) 

Also, we may use the rectangle lattice of scales and one-dimentional wavelet decomposition : 

f(x u x 2 )= < /.*i,/®*j,fc > *j,/®*j,k(&i,&2), (57) 

i,t;j,k 

where bases functions 

*M <8> *j,fc (58) 

depend on two scales 2~ l and 2~- 7 . 

After construction the multidimensional bases we obtain our multiscale/multiresolution representations for 
observables (operators, symbols), states, partitions via the variational approaches [2]- [8] as for c-BBGKY as for its 
quantum counterpart and related reductions but before we need to construct reasonable multiscale decomposition 
for all operators included in the set-up. 

2.3. FWT Decomposition for Observables 

One of the key point of wavelet analysis approach, the so called Fast Wavelet Transform (FWT) , demonstrates 
that for a large class of operators the wavelet functions are good approximation for true eigenvectors; and 
the corresponding matrices are almost diagonal. FWT gives the maximum sparse form for wide classes of 
operators [10]. So, let us denote our (integral/differential) operator from equations under consideration as T 
{L 2 {R n ) — ► L 2 {R n )) and its kernel as K. We have the following representation: 

<Tf l9 >=J f K(x,y)f(y)g(x)dxdy. (59) 

In case when / and g are wavelets ipj t k = 2^ 2 ip(2^x — fc), (21) provides the standard representation for operator 
T. Let us consider multiresolution representation . . . C V 2 C V\ C Vo C V-\ C V- 2 .... The basis in each Vj 
is (fij t k(x), where indices fc, j represent translations and scaling respectively. Let Pj : L 2 (R n ) — > Vj (j G Z) be 
projection operators on the subspace Vj corresponding to level j of resolution: (Pjf)(x) — J^k < /)Vj,fc > ¥j,k(x). 
Let Qj = Pj-i — Pj be the projection operator on the subspace Dj (Vj-i = Vj (B Dj), then we have the following 
representation of operator T which takes into account contributions from each level of resolution from different 
scales starting with the coarsest and ending to the finest scales [10]: 

T = zZiQiTQi + Qi TP 3 + PjTQj). (60) 

jez 

We need to remember that this is a result of presence of affine group inside this construction. The non-standard 
form of operator representation is a representation of operator T as a chain of triples T — {Aj , Bj , Tj }jez 1 acting 
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on the subspaces Vj and Df Aj : Dj — > Dj,Bj : Vj — ► Dj,Tj : Dj — ► V}, where operators {Aj, Bj,Tj}j € z are 
defined as = QjTQj, Bj = QjTPj, Tj = PjTQj. The operator T admits a recursive definition via 



Tj = 7 + ^ , (61) 



where Tj — PjTPj and Tj acts on Vj : Vj — > Vj. So, it is possible to provide the following "sparse" action of 
operator Tj on elements / of functional realization of our space of states H: 

(Tjf)(x) = ]T [2-'" nfj,k-e) VjM (62) 

in the wavelet basis ipj. k {x) = 2~i/ 2 ip(2~i x — k), where 

fi,k-i = 2^' /2 / f(xM2~ j x -k + £)dx (63) 



are wavelet coefficients and re are the roots of some additional linear system of equations related to the "type 
of localization" [10]. So, we have the simple linear parametrization of matrix representation of our operators in 
localized wavelet bases and of the action of this operator on arbitrary vector/state in proper functional space. 

2.4. Variational Approach 

Now, after preliminary work with (functional) spaces, states and operators, we may apply our variational ap- 
proach from [2]- [8]. 

Let L be an arbitrary (non)lincar differential/integral operator with matrix dimension d (finite or infi- 
nite), which acts on some set of functions from L 2 (f2®"): \& = ty(t, Xi, X2, ■ ■ •) = (^! 1 (t,Xi,X2, ■■■),■■ 

Xi,X2, •••))' ii 6 O C R 6 , n is the number of particles: 

T* = L(Q,t,Xi)V(t,Xi) = 0, (64) 
Q = Qd ,di,d2,...(t,x 1 ,x 2 , ■ ■ ■ ,d/dt,d/dx 1 ,d/dx 2 , ■ ■ ■ ,J (in) 

do ,di ,d,2 , - - ■ Q io Q i\ Q i"2 f 

j0iUi»2,-.. = l 

Let us consider now the N mode approximation for the solution as the following ansatz: 

N 

V N (t,xi,X2,...)= ^2 a io i 1 i 2 ..A io <g>Bi 1 <g,C i2 ...(t,xi,X2,...) (65) 

•Oi'li*2r"-l 

We shall determine the expansion coefficients from the following conditions (related to proper choosing of vari- 
ational approach): 

Cki,k 2 ,... = J(L* N )A ko (t)B kl (x 1 )C k2 (x 2 )dtdx 1 dx 2 ... = (66) 

Thus, we have exactly dN n algebraical equations for dN n unknowns a,i 0i i lt .... This variational approach reduces 
the initial problem to the problem of solution of functional equations at the first stage and some algebraical 
problems at the second one. It allows to unify the multiresolution expansion with variational construction [2] -[8]. 

As a result, the solution is parametrized by the solutions of two sets of reduced algebraical problems, one 
is linear or nonlinear (depending on the structure of the generic operator L) and the rest are linear problems 
related to the computation of the coefficients of reduced algebraic equations. It is also related to the choice of 
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exact measure of localization (including class of smothncss) which arc proper for our set-up. These coefficients 
can be found by some functional/algebraic methods by using the compactly supported wavelet basis functions 
or any other wavelet families [10]. 

As a result the solution of the equations/hierarchies from Section 1, as in c- as in q- region, has the following 
multiscalc or multiresolution decomposition via nonlinear high-localized eigenmodes 

W(t,xi,x 2 ,...)= o-ijU 1 ®V 3 {t lXu . . 

yj{t) = v j,slo W{t) + yi^i), ^ ~ 2 l , (67) 

1>N 

U i (x s ) = Uti low (x s )+ J2 UL(k s m x s ), k' m ~2 m , 

m>M 

which corresponds to the full multiresolution expansion in all underlying time/space scales. The formulae (67) 
give the expansion into a slow part and fast oscillating parts for arbitrary N, M. So, we may move from the 
coarse scales of resolution to the finest ones for obtaining more detailed information about the dynamical process. 
In this way one obtains contributions to the full solution from each scale of resolution or each time/space scale or 
from each nonlinear eigenmode. It should be noted that such representations give the best possible localization 
properties in the corresponding (phase)space/time coordinates. Formulae (67) do not use perturbation techniques 
or linearization procedures. Numerical calculations are based on compactly supported wavelets and wavelet 
packets and on evaluation of the accuracy on the level N of the corresponding cut-off of the full system regarding 
norm (45): 

\\W N+1 ~W N \\<e. (68) 
3. MODELING OF PATTERNS 

To summarize, the key points are: 

1. The ansatz-oriented choice of the (multidimensional) bases related to some polynomial tensor algebra. 

2. The choice of proper variational principle. A few projection/ Galerkin-like principles for constructing 
(weak) solutions are considered. The advantages of formulations related to biorthogonal (wavelet) decomposition 
should be noted. 

3. The choice of bases functions in the scale spaces Dj from wavelet zoo. They correspond to high-localized 
(nonlinear) oscillations/excitations, nontrivial local (stable) distributions/fluctuations, etc. Besides fast conver- 
gence properties it should be noted minimal complexity of all underlying calculations, especially in case of choice 
of wavelet packets which minimize Shannon entropy. 

4. Operator representations providing maximum sparse representations for arbitrary (pseudo) differential/ 
integral operators df/dx, d n f/dx n , J T(x,y)f(y)dy), etc. 

5. (Multi)linearization. Besides the variation approach we can consider also a different method to deal with 
(polynomial) nonlinearities: para-products-likc decompositions. 

To classify the qualitative behaviour we apply standard methods from general control theory or really use 
the control. We will start from a priori unknown coefficients, the exact values of which will subsequently be 
recovered. Roughly speaking, we will fix only class of nonlinearity (polynomial in our case) which covers a 
broad variety of examples of possible truncation of the systems. As a simple model we choose band-triangular 
non-sparse matrices (a^ ). These matrices provide tensor structure of bases in (extended) phase space and 
are generated by the roots of the reduced variational (Galerkin-like) systems. As a second step we need to 
restore the coefficients from these matrices by which we may classify the types of behaviour. We start with the 
localized mode, which is a base mode/eigenfunction, (Fig. 1, 9, 10 from Part I), corresponding to definitions 
from Section 2.2, Part I, which was constructed as a tensor product of the two Daubechies functions. Fig. 1, 4 
below demonstrate the result of summation of series (67) up to value of the dilation/scale parameter equal to 
four and six, respectively. It's done in the bases of symmlets [10] with the corresponding matrix elements equal 
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to one. The size of matrix of "Fourier-wavelet coefficients" is 512x512. So, different possible distributions of the 
root values of the generical algebraical systems (66) provide qualitatively different types of behaviour. Generic 
algebraic system (66), Generalized Dispersion Relation (GDR), provide the possibility for algebraic control. The 
above choice provides us by a distribution with chaotic-like equidistribution. But, if we consider a band-like 
structure of matrix (a^ ) with the band along the main diagonal with finite size (<C 512) and values, e.g. five, 
while the other values are equal to one, we obtain localization in a fixed finite area of the full phase space, i.e. 
almost all energy of the system is concentrated in this small volume. This corresponds to waveleton states and is 
shown in Fig. 2, constructed by means of Daubechies-based wavelet packets. Depending on the type of solution, 
such localization may be conserved during the whole time evolution (asymptotically-stable) or up to the needed 
value from the whole time scale (e.g. enough for plasma fusion/confinement in the case of fusion modeling by 
means of c-BBGKY hierarchy for dynamics of partitions). 

4. CONCLUSION 

So, by using wavelet bases with their best (phase) space/time localization properties we can describe the localized 
(coherent) structures in quantum systems with complicated behaviour (Fig. 1,4). The modeling demonstrates 
the formation of different (stable) pattern or orbits generated by internal hidden symmetry from high-localized 
structures. Our (nonlinear) eigenmodes are more realistic for the modelling of nonlinear classical/quantum dy- 
namical process than the corresponding linear gaussian-like coherent states. Here we mention only the best 
convergence properties of the expansions based on wavelet packets, which realize the minimal Shannon entropy 
property and the exponential control of convergence of expansions like (67) based on the norm (45). Fig. 2 cor- 
responds to (possible) result of superselection (einselection) [1] after decoherence process started from entangled 
state (Fig. 5); Fig. 3 and Fig. 6 demonstrate the steps of multiscale resolution during modeling of entangled 
states leading to the growth of degree of entanglement. It should be noted that we can control the type of 
behaviour on the level of the reduced algebraical variational system, GDR (66). 
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Figure 1. Level 4 MR A. 



Figure 4. Level 6 MRA. 
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Figure 2. Localized pattern, (waveleton) Wigner function. 



Figure 5. Entangled- like Wigner function. 
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Figure 3. Interference picture on the level 4 approximation Figure 6. Interference picture on the level 6 approximation 
for Wigner function. for Wigner function. 
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